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We present a complete next to next to leading order QCD corrections to the resonant 
sneutrino and charged slepton production cross sections at hadronic colliders such 
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violating super symmetric model. We have demonstrated the role of these corrections 
in reducing uncertainties resulting from renormalisation and factorisation scales and 
thereby making our predictions reliable. We have incorporated soft gluon effects at 
N^LO level in order to study the stability of our results under perturbation. The 
NNLO corrections are found to be large and significant. The results obtained in 
this article are also applicable to resonance production of any color-neutral scalar. 
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1 Introduction 



The Standard Model (SM) of particle physics is undoubtedly a very successful model even 
though its scalar sector remains in an unsatisfactory state of affairs. This is because of the 
fact that the Higgs boson in the scalar sector which is responsible for symmetry breaking 
mechanism is still missing. There are also other issues such as gauge hierarchy problem, 
dark matter, baryogenesis, gauge unification etc., that are still not fully understood within 
the framework of the SM, pointing to physics beyond the standard model. In the past, 
the above mentioned issues were being addressed by going beyond the SM. Two of the 
most attractive classes of such models contain features incorporating supersymmetry [1] 
and/or grand unification [2] (scenarios with a low intermediate scale [3]). These models 
contain large number of new particles including new scalars. The couplings of the first 
generation SM fermions to these new scalars need not be suppressed. This provides us a 
new platform to study their potential new physics signals. 

The electroweak gauge invariance ensures both baryon {B) and lepton (L) number 
conservations within the SM, but this is not the case for supersymmetry (SUSY). The 
most general superpotential respecting the gauge symmetry of the SM contains bilinear 
and trilinear terms which do not respect B ot L conservations. A discrete symmetry 
called i? -parity can be used to forbid such terms. The corresponding conserved quantum 
number is given by Rp = (_i ^^b+l+ss^ where S is the spin of the particle. In the minimal 
supersymmetric standard model (MSSM), this symmetry has been originally imposed to 
suppress the rapid proton decay, it also guarantees stability of lightest supersymmetric 
particle (LSP) and hence a natural candidate for cold dark matter. 

The introduction of this symmetry is not the only option to suppress proton decay, 
there exists several alternatives through introduction of other symmetries. Hence, it is 
of phenomenological interest to consider possible violations of R-parity and study their 
experimental consequences. 

The possible i?-parity violating ( Jlp) terms in the superpotential can be parametrised 

as 

W^, = fi^UH, + A,,,L,L,i?^ + K,,UQ,Dl + K^kU^Dpi , (1) 

where Li and Qi are the SU (2)-doublet lepton and quark superfields, E'[, f/f , the sin- 
glet superfields and Hi the Higgs superfields. Note that Xijk is antisymmetric under the 
interchange of the first two indices and X"j^. is antisymmetric under the interchange of 
the last two. The first three terms in eqn.(l) violate L and the last term violates B 
conservation. We need to have at least one of the two sets of couplings to be vanishingly 
small in order to satisfy the constraints coming from the non-observance of proton decay. 
In the following, we assume that i? is a good symmetry of the theory, or in other words 
all of X"j}^ are zero. This can also suppress all dimension six operators leading to proton 
decay (see [4]) along with the dimension five ones. In [4,5], this scenario has been moti- 
vated within certain theoretical frameworks and makes simpler, problem of preservation 
of GUT-scale baryon asymmetry [6]. The presence of the other terms can affect the 
baryon asymmetry of the universe. Fortunately these bounds are highly mo del- dependent 
and hence we evade them (see [7]). 
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Each interaction in eqn.(l) has its unique experimental signatures through low-energy 
phenomenology and/or in resonance production at colliders. Large number of bounds on 
these coupling result from the studies of low and intermediate energy processes, namely 
rare decays with lepton and/or hadron flavor violation. The neutral and charged current 
universality in lepton and quark sectors, high precision measurements in anomalous mag- 
netic and electric dipole moments and observables with CP violation provide stringent 
constraints on these couplings at low energies. Terms involving Xij^ lead to resonant sncu- 
trino production in e+e" collider [8, 9] and those involving X'-jj. lead to resonant squark 
production in hadron-hadron collisions [10,11]. Similarly, the X'^^jj^ terms can contribute 
to both resonant squark production at an e^p coUiders [12] as well as to resonant charged 
slepton and sneutrino production at hadron colliders [10,13]. A simplified search strategy 
at colliders could be to assume the existence of one dominant R-parity violating coupling 
at a time. In the past most of the studies rely on this assumption. Since the absence 
of tree-level flavor changing neutral current processes lead to severe constraints on the 
simultaneous presence of more than one X[j,^ [14], we shall restrict ourselves to only one 
non-zero A^^j.. In the rest of the paper, we will concentrate only on the X\jf, couplings that 
can affect the resonant production. That part of the Lagrangian can be written in terms 
of the component fields as 



The squarks in supersymmetric theories behave as leptoquarks and the charged slep- 
tons/sneutrinos behave as charged/neutral Higgses of multi-Higgs-doublct scenario. Hence, 
one would expect non-zero values for these couplings will have important phenomenolog- 
ical consequences. 

Quantum Chromo dynamics (QCD) plays an important role in hadron colliders as the 
underlying scattering processes involve quarks (anti-quarks) and gluons. The leading 
order (LO) scattering processes involving quarks (anti-quarks) and gluons often predict 
results that are sensitive to large theoretical uncertainties through non-perturbative par- 
ton densities and missing higher order perturbative radiative corrections. In the next 
sections, wc will elaborate on the sources of these uncertainties and provide systematic 
methods to reduce them. In [15], first results on the next-to-leading order (NLO) QCD 
corrections to sneutrino and charged slepton productions at hadron colliders were re- 
ported. These results were later confirmed and further NLO SUSY-QCD corrections were 
systematically incorporated in the analysis by the authors of [16]. It was found that 
the NLO QCD effects were quite large ~ 10% — 40% at both Tevatron as well as LHC 
and hence they were used by both CDF [17] and DO [18] collaboration to analyse their 
data (Run-1 as well as Run-11 data). In their analysis to set bound on these R-parity 
violating couplings, K-factor for SM background [] was considered at the next to next to 
leading order (NNLO) level while for the R-parity violating effects only NLO K-factor 
was used. Therefore, it is desirable to compute the K-factors for the resonant sneutrino 
and/or charged slepton productions at NNLO in QCD. These results will quantitatively 
improve the analysis based on high statistics data available in the ongoing and future ex- 
periments. Prom the theoretical point of view, higher order radiative corrections provide 




(2) 



dkRiiLUjL - dkRUjdiL - {.(-ihYUjLdlR + h.c. 
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a test of the convergence of the perturbation theory and hence the rehable comparison 
of data with the theory predictions is possible. We have also studied the effects of soft 
gluons which often dominate at hadronic collisions. This also opens up the possibility of 
resumming them to all orders in perturbation theory through suitable framework. The 
fixed order perturbative results most often suffer from large uncertainties due to the pres- 
ence of renormalisation and factorisation scales. They get reduced as we include more and 
more terms in the perturbative expansion thanks to renormalisation group invariance. It 
is also important to estimate the uncertainties coming from the choice of parton density 
sets available in the literature. Due to these reasons, we think that it is worthwhile to 
perform the next-next-to-leading order (NNLO) calculation to the aforementioned pro- 
cesses. In our previous paper [15], we used fixed ^pCOuphng (A' = 0.01) based on the 
assumption that its scale dependence is very weak, in addition the value of the coupling 
A' was immaterial and only served to set an overall normalisation for the cross-section. 
In the present paper we have systematically included its scale dependence through the 
renormalisation group equations and we will discuss the impact of it in the next sections. 

We briefly review the constraints on A' from low-energy phenomenology. Non-zero 
A's can lead to additional four-fermion operators that may contribute to meson decays, 
neutral meson mixing, some of which may be forbidden otherwise. In Table 1, we list the 
currently known bounds on several of these couplings^. The strongest bound is on A^^^^ and 
is derived from non-observation of neutrinoless double beta decay (a) [20]. The others are 
much weaker and are derived from (6) upper bound on the mass of the i^e [5,21-23]; data 
on (c) charged- current universality [9]; (d) atomic parity violation [24]; (e) r ■nvr and 
D — >■ Klv [22]; and (/) D^-D^ mixing [14]. Since these bounds are derived from effective 
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0.001 


211 


0.09^^ 


311 


0.10") 


112 


0.02 ^) 


212 


0.09 ^) 


312 


0.10") 


121 


0.035 


221 


0.18") 


321 


0.20 ^) 


122 


0.02 


222 


0.18") 


322 


0.20 ^'^ 



Table 1: The upper bounds on the A'-type couplings of interest for a common sfermion 
mass m = 100 GeV. The superscripts refer to the specific experiments leading to the 
constraints and as described in the text. 

four-fermion operators, they typically scale like the mass of the exchanged sfermion^. 
These bounds are actually applicable only to particular combinations of couplings and 
masses and reduce to those in the table only under the assumption of only one coupling 
being non-zero. In addition, in meson decays, most often it is the squark that is exchanged; 
hence charged sleptons/sneutrinos could very well be much lighter without contradicting 
the bounds. 

"'A more complete list can be found in refs. [19]. 

^Of those listed in Table 1, the only exceptions to this rule are the bounds for A'^ and A'122 [5,20-22]. 
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In the next section, we describe the computation of resonant production of scalar/pseudo 
scalar to NNLO in perturbative QCD. The results are obtained in the MS scheme. We 
then proceed to study the impact of our results on their production cross sections for both 
Tevatron and LHC energies. Finally we summarise our findings in the conclusion. Our 
analytical results that go into the numerical code are presented in the Appendix. 



2 Computation of partonic coefficient functions to 
order ai 



In this section, we describe in detail, the computation of second order (ct^) QCD radiative 
corrections to resonant production, in hadron colliders, of a neutral scalar particle 0(^) 
which couple to fermionic fields through Yukawa interaction given by the action: 

s = I A„,(o = x'l d'^ m mm (3) 

where A' is the coupling strength of the interaction. We present our results in such a 
way that they can be used for a similar study on charged scalar production as well and 
hence are applicable to a detailed study on resonant production of sneutrinos and charged 
sleptons which is the main goal of this present work. The inclusive hadronic cross section 
for the reaction 



H,{P,) + H2{P2) ^ (1>{P5) + X , (4) 



is given by 



4> (i^r) dxi dx2 , 2 x . / 2 x a ^ ^ 222 

^^-^ a,b=q,q,g ^1 "^-^/^l ^2 \X1X2 

2 

TV) 

with r = -^ , S^{P, + P2f , Pl^ml, (5) 

where Hi and H2 denote the incoming hadrons and X represents an inclusive hadronic 
state. The parton densities denoted by fc{xi,fip) (c = q,q,g) depend on the scaling 
variables Xi {i — 1,2) through pi — XiPi and the mass factorization scale /ip- Here pi {i — 
1, 2) are the momenta of incoming partons namely quarks, anti-quarks and gluons. The 
coupling constant A' gets renormalised at the renormalisation scale fifi due to ultraviolet 
singularities present in the theory. The factorisation scale is introduced on the right hand 
side of the above equation to separate long distant dynamics from the perturbatively 
calculable short distant partonic coefficient functions Aab- ^ab depends on both //^ and 
in such a way that the entire scale dependence goes away to all orders in perturbation 
theory when convoluted with appropriate parton densities. This is due to the fact that 
the observable on the left hand side of the above equation is renormalisation group (RG) 
invariant with respect to both the scales. This implies 
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0. 



(7) 



The partonic coefficient functions that appear in eqn.(5) are computable in perturbative 
QCD in terms of strong coupUng constant gg- The ultraviolet singularities present in 
the theory are regularised in dimensional regular isat ion and are removed in MS scheme, 
introducing the renormalisation scale at every order in perturbative expansion. In 
addition, the Yukawa coupling A' also gets renormalised due to strong interaction dynam- 
ics. Hence, for our computation, we require only two renormalisation constants to obtain 
UV finite partonic coefficient functions, Aab- These constants are denoted by Z{fifi) and 
Zy{f^R), where the former renormalises the strong coupling constant gs and the later 
Yukawa coupling A'. 

We define the bare strong coupling constant by as = ^f/lGvr^, gs being the dimen- 
sionless strong coupling constant in n = 4 + e, with n being the number of space time 
dimensions. The bare coupling constant cig is related to the renormalised one, as{n'ji) 
the following relation: 



Setts = Z{as, /x^ Hr, €)asinR) 



^R. 



2 \ 2 



(8) 



The scale /i comes from the dimensional regularisation in order to make the bare coupling 
constant cjs dimensionless in n dimensions. S^ is the spherical factor characteristic of 
n- dimensional regularisation. 

The renormalisation constant that relates the bare coupling constant ds to the renor- 
malised one as{iJ,%) through the eqn.(8) is given by 



Z(as,//^//^,e) = 1 + ds 



s. 



2/3o 



-|- a 
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+ o{dt) . 



(9) 



The coefficients /3j for i = 1, 4 can be found in [25] for 5*^7 (N) QCD expressed in terms 
of the color factors: 



2N ' ' 2 

and n/ is the number of active flavors. Similarly, for the Yukawa coupling, we have 



(10) 



SeX' = Zx'{ds, /i^ fiR, e)A'(/i^) 



2X1 



(11) 



where 



Zy{ds,iJi ,ijLji,e) = 1-FaJ— Se 



a- 



+ di 



2 Po 70 



5 



+ -171 



ifi(7c)=-4A(7o)%|/J„S; 



1 / 2 ^ 8 ^ \ 1 /2 

^ 2 70 71 - 3 ^70-3^0 71 3 



+ 0{at) . (12) 



The anomalous dimensions ji for i = 1, 4 can be obtained from the quark mass anoma- 
lous dimensions given in [26]. While the C(a^) terms in both Z and Zy do not contribute 
to partonic sub processes computed to order a^, they determine the scale evolution of 
both tts and A' to NNLO through renormalisation group equations: 

d °° 

"/^R i=l 

J oo 

/^^^lnA'(/.^) = -Ea:(/.|)7.-i. (13) 

and constitute dominant soft gluon contribution to order a^. The perturbatively calcula- 
ble Aab can be expanded in powers of strong coupling constant as(A*ie) ^ 

oo 

Aa6 (x, ml, III, fil) = J2 O^'sitAd^^'al {^^ ^% I^F, lAi) 

i=0 

Aab gets contributions from various partonic reactions. 

Having studied the UV renormalisation constants relevant for our computation, we 
now hst out various partonic sub processes that will contribute to NNLO order (a^). The 
leading order partonic reaction (fig.(l)) is given by 

qi + qj^(l) (14) 

and its contribution is found to be proportional to 5(1 — a;) where x — ml/s with partonic 
center of mass energy s — {pi +^2)^- At NLO (fig. (2, 3,4)), we have 

Qi + Qj + one loop, 
Qi + Qj (p + 9, 

<li{li)+9 (t> + (lMj)- (15) 

The subscripts i,i in the quark (anti-quark) fields denote the flavor indices. Hi = j, then 
(p will be neutral scalar otherwise it will denote a charged scalar. 

At NNLO, there are several new channels open up and we discuss them one by one: 

• processes with quark anti-quark pair in the initial state with no gluon (fig. (5)), one 
gluon (fig. (6)) and two gluons (fig. (7)) in the final state along with 0: 

Qi + Qj + two loop , 
Qi + Qj + 9 + one loop , 

Qi + q^ ^ (f) + 9 + 9- (16) 
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• processes (fig. (8)) with quark anti-quark pair in the final states along with 0: 

qi + Qj ^ (p + qk + Qi, (17) 
and quark (anti-quark) quark (anti-quark) pair in the final state with (p (fig. (9)): 

qi + Qj ^ (p + Qk + qi, 

qi + qj ^ (p + qk + qi- (is) 

The contributions coming from the last two reactions (eqn.(18)) will have two pos- 
sibilities: the final states with identical quarks (anti-quarks) or non-identical quarks 
(anti-quarks). For identical quarks (anti-quarks) in the final state (fig. (10)), the 
number of t channel processes doubles as we need to include processes with final 
state quarks (anti-quarks) interchanged. We also need to appropriately multiply the 
statistical factor 1/2. 

• processes with quark (anti-quark) and gluon in the initial states 

qi{qi) + 9 (I) + qj{qj)+ one loop, 

qiiqi) + 9 (l) + qAqj) + 9- (i9) 

• processes with pair of gluons in the initial state 

g + g^^ + qi + q^. (20) 



The calculation of various contributions from the partonic reactions involves careful 
handling of divergences that result from one [27] and two loop [28] integrations in the 
virtual processes and two and three body phase space integrations in the real emission 
processes. The loop integrals often give ultraviolet, soft and collinear divergences. But 
the phase space integrals give only soft and collinear singularities. Soft divergences arise 
when the momenta of the gluons become zero while the collinear diverges arise due to 
the presence of massless partons. We have regulated all the integrals in dimensional 
regularisation with space time dimension n — 'i + e. The singularities manifest themselves 
as poles in e. 

We have reduced all the one loop tcnsorial integrals to scalar integrals using the 
method of Passarino-Veltman [29] in 4 + e dimensions and evaluated resultant scalar 
integrals exactly. The form factor that contributes to reactions in eqn (16) is obtained 
using the dispersion technique [30] and is presented in the appendix. Two and three 
body phase space integrals are done by choosing appropriate Lorentz frames [31]. Since 
we integrate over the total phase space the integrals are Lorentz invariant and therefore 
frame independent. Several routines are made using the algebraic manipulation program 
FORM [32] in order to perform tensorial reduction of one loop integrals and two and three 
body phase space integrals. 

The UV singularities go away after performing rcnormalisation through the constants 
Z and Zy. The soft singularities cancel among virtual and real emission processes [33] at 



7 



every order in perturbation theory. The remaining coUinear singularities are renormahsed 

systematically using mass factorisation [34] as follows. For more details on the compu- 
tation of NNLO QCD corrections to process of the kind considered here can be found 
in [35]. Let us denote the resulting UV and soft finite partonic cross sections as 

A„6 (x, ml, /i^) = x'\iJ,'^) ^"'''^ ^^-^ ' 



where 



m=3 i=3,i^5 v / 
m 

3=1 

Mab^X <t> = {pi,P2\OiO)\X,p5) with \X,P5) = \p3,P4,Pq- ■ ■pm,P5) , 

(21) 

0(0) is the interaction operator responsible for the reactions with scalars and Kab repre- 
sents the spin and colour average over the initial states. Mass factorisation allows us to 
express the coUinear singular partonic cross section Aab in terms of pair of singular tran- 
sition functions Tcdix, e), namely Altarelli-Parisi kernels and finite partonic coefficient 
function Aa^: 



c,d=q,q,g "^^ ^1 "^^'^i ^2 

xA,J^,ml,fil,fil) . (22) 
The transition functions are perturbatively calculable in powers of as(/x|'): 



oo 



r„.(a:,/.^,6) = E«:(/^F)rlV(^,e) ■ (23) 

i=0 



In MS mass factorisation scheme, they are found to be (suppressing the arguments x and 

ri? = 5„,5(i-x), (24) 

^'2 - (25) 



and 



= 2^ E (^i? ® + 2/3oPi?) + ^i^i^^ ■ (26) 



8 



More explicitly we have: 



T^{2),NS 



qq 



r(2),5/- 



.{2),NS 
QQ 



p(2),JVS _ j_p(l)- 

^99 — 2e^^ 



p(2),5/- _ p(2),5/- _ p(2),5/- 
J- gg — ^ qq ^ ^ QQ 



— p(°) (5^) p(o) + — p(i)'-^/- 



r(2) ^ 



99 



99 



2^2 1 99 ^ 99 ^ 99 ^ 99 ^ ^U-' qg ( ^ 2^ Q9 



In order to determine Aab, we set /jr = fip and expand Aa6 (similarly Aab) as 



\ /^F/ i=0 \ A'-F/ 



(27) 



(28) 



Substituting the above equation (eqn. (28)) in the eqn.(22), we obtain to order as{iip) 





= A^'J 


QQ 


QQ 


A(i) 

Q9 


= A« 

99 



AW--Pio)®Ai°J 



99 



99 



A(i)-ipW®A(°J 

99 ^ 99 ^ 99 



(29) 



and to order aKiij^) 



99 



- 2r« ® a(°) 



99 
1 



99 



-PW® Aiy--Pio)®A« 



99 



99 



99 



99 



p(0) ^ p(0) ^ ^(0) 



c2 99 



99 



99 



A(2) 

99 



c2 QQ 



A(2) 

99 



A(2) - lp(o) <^ aW - -Pjo) (g) PW (g) Ai°^ 



99 



99 



99 



,,2 99 



99 



99 



(30) 



We have computed Aab for all those reactions that contribute to scalar production to order 
and substituted them into eqns. (29,30) to obtain finite partonic coefficient functions 
Aab- The results are lengthy and hence they are presented in the Appendix after setting 
l-^R = f^F- It is straightforward to obtain the fi^ dependence of the partonic coefficient 
functions presented in the Appendix using the RG equations given in eqns. (7). The mass 
factorisation implies DGLAP evolution equation: 



^^F-^fa{x, /^f) = I X] ^cib {x, (8) fb{x, n%) , 



(31) 
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where Pah are Altarelli-Parisi splitting functions [36]. The DGLAP evolution equation 

determines the scale evolution of the parton densities appearing in eqn.(5). Note that the 
RG equations of strong, Yukawa couplings (eqn.(13)) along with the above DGLAP evolu- 
tion equations (eqn.(31)) control both factorisation and renormalisation scale dependences 
of at every order in perturbation theory. 

We have made several checks on our NNLO results, both analytically as well as nu- 
merically. First and the foremost check is the observation of cancellation of all the poles 
in e that result from UV, soft and coUinear divergences after all the appropriate renor- 
malisation constants and factorisation kernels are systematically taken into account. The 
second check involves the comparison of our results against those computed for Higgs 
production through bottom quark annihilation [37]. Note that we have presented our 
results in such a way that they can be used for both neutral scalar (sneutrino) and 
charged scalar (sleptons) productions. So, the comparison against [37] is possible only 
after combining various pieces of the coefficient functions and then multiplying x and 
setting Ca = "i-iCp — 4:/3,Tf — 1/2 in our expressions. We found complete agreement 
with [37] which serves as an important check on our computation. Finally, we have re- 
produced all those plots given in [37] using our code in order to check the correctness of 
numerical code. In the next section we will discuss in detail the numerical impact of our 
results for both sneutrino and charged slepton productions at hadron colliders. 

Recently, there have been several breakthroughs in understanding the structure of 
perturbative series to all orders in perturbation theory, thanks to explicit results on form 
factors and Altarelli-Parisi sphtting functions (see [38-43]) to three loop level in QCD. The 
computation of such quantities reveals the long distance physics resulting from the soft 
gluon emissions in the scattering processes. The soft gluon contributions to hadronic cross 
sections often dominate over the rest in the region where the partonic scaling variable x 
approaches unity. Due to the peculiar behavior of the partonic fluxes in this region, often 
these effects need to be resummed to all orders in perturbation theory. Resummation of 
soft gluons for hadronic reactions can now be achieved upto next to next to next to leading 
logarithm (N^LL) level using the available three loop results (see [44-47]). More on the 
structure of perturbative results both in fixed order as well as in the resummed quantities 
can be found in [47-50] . In [47] , soft gluon enhanced partonic contributions were obtained 
for Drell-Yan, Higgs productions at hadronic coUiders. This was achieved by using the 
coUinear factorisation property of scattering cross sections, the Sudakov resummation of 
soft gluon effects and applying various renormalisation group techniques. In [51], order 
soft gluon contributions to Higgs production through bottom quark annihilation process 
at hadron colliders were presented for the first time using the soft gluon enhanced cross 
sections thus obtained. Since the coupling of SM Higgs boson to bottom quarks and 
that of sneutrino/sleptons are both of Yukawa type, we can use of [51] to study 
the soft gluon effects on the snuetrino and charged slepton production cross sections at 
hadron colliders. We will present the numerical impact of these effects towards the end of 
next section. Such a study presents a quantitative estimate on the missing higher order 
contributions to the processes under study. 
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3 Results and Discussion 



Having obtained the compact analytic results for the partonic coefficient functions /S.ab 
for various subprocesses to NNLO in perturbative QCD, we now endeavor to study their 
impacts on the resonant production of sneutrino and charged slepton at the LHC {s/S = 
14 TeV) and for the Run II of Tevatron (v^ = 1.96 TeV). As discussed in the Introduction, 
we will limit ourselves only to contributions from the first generation of quarks . Since 
at hadron colliders, the resonant production is through the interaction term X'^jj^LiQjD^ 
in the Lagrangian (see eq.(l)), for j, k = 2,3, the production rate will be suppressed due 
to the low flux of the sea quarks. To obtain the production cross section to a particular 
order, one has to convolute the partonic coefficient functions Aab with the corresponding 
parton densities /„, both to the same order. Further the coupling constants as{iJ,R) and 
X'{IjLr) should also be evaluated using the corresponding RGEs (eqn.(7)) computed to the 
same order. 

Following the Ref [26], the solution to RGE (second equation in eq.(7)) for A'(/x|j) is 
given by , 



with 



C{as) = af^Y.<C,. (33) 

1=0 

The Ci are given by 

Co = 1, Ci = Ai, 

C2 = ^{Al + A2), Cs = ^{Al + 3A,A, + 2A,), (34) 

with 

Aq = Co, Ai = ci- 6iCo, ^2 = C2 - 6iCi + co{bl - 62), 

A3 = C3-biC2 + Ci{bl-b2)+Co{bib2-bi{bl-b2)-b3), (35) 



and 



Ci = bi^ (36) 



and /^o is some reference scale at which both as well as A' are known. We have 
numerically evaluated as(/x|j) and A'(/x|j) to relevant order namely LO, NLO and NNLO 
by truncating the terms in the RHS of eqn.(13). We have used A'(M|) = 0.01 irrespective 
of flavor and mass of the sncutrino/charged slepton. We have used the latest MSTW 
parton densities [52] in our numerical code and the corresponding values of as{Mz) for 
LO, NLO and NNLO provided with the sets. 
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3.1 Sneutrino Production 



The total sneutrino production cross section as function of its mass is plotted in fig. 11 for 
LHC (left panel) and Run 11 of Tevatron (right panel) energies. We have set the renormal- 
isation scale to be the mass of the sneutrino, fiR = rric,. The pair of lines corresponds to 
the two extreme choices of factorisation scale: iip = iTii> (upper) and iip = O.lmp (lower). 
The plots clearly demonstrate that the NNLO contributions reduce the factorisation scale 
dependence improving the theoretical predictions for sneutrino production cross section. 
To quantify the percentage variation with respect to the factorisation scale we define 

Sa'^iUp) = = /^fl = iTT'o) - (/^F = O.lmp, /ir = nip) ^^^^ 

ilJ'F ^ IJ'R^ mo) 

where / = LO, NLO, NNLO. The cross section falls off with the sneutrino mass due to 





Mass Range (m/>) 


LO in % 


NLO in % 


NNLO in % 


Tevatron 


100 GeV to 1 TeV 


-4.6 to 140 


-2.2 to -39 


-1.1 to -10.6 


LHC 


100 GeV to 450 GeV 
510 GeV to 1 TeV 


40 to 2.3 
-1.2 to -17 


23 to 1.8 
-0.22 to -5.2 


9.3 to 1.1 
-0.24 to -1.2 



Table 2: da^'^/ip) in % in the mass range between 100 GeV and 1 TeV at Tevatron and 
LHC energies. 



the availability of phase space with respect to the mass, the choice of fi^ = mp and the 
parton densities. The latter effect, understandably, is more pronounced at the Tevatron 
than at the LHC. 

In order to estimate the magnitude of the QCD corrections at NLO and NNLO, we 
define the K-factors as follows: 

K^'^ = 4lo/4o = 4nlo/4o- 

In fig. 12, we have plotted both K^'^^ (i = 1, 2) as a function of sneuttrino mass. We have 
chosen jjip — n^ — mo for this study. At the LHC, The K^^'> varies between 1.23 to 1.46 
and between 1.27 to 1.52 in the mass range 100 GeV < m < 1 TeV . At the Tevatron, 
we find that K^^"* varies between 1.55 to 1.53 and fT^^^ between 1.65 to 1.85 for the same 
mass range. Note that numbers for K^^"* differ slightly from those given in our earlier 
work [15] due to the running of A' in the present analysis. We also observe that K factor 
is much bigger at the Tevatron compared to that of at the LHC. The reason behind this 
is attributed to the different behavior of parton densities at the Tevatron and the LHC. 
Note that parton densities rise steeply as a; — )■ and fall off very fast 1, which 

means the dominant contribution to the production results from the phase space region 
where x ~ t{— ml/S) becomes small, r at Tevatron (0.05 ^ r ^ 0.5) is larger compared 
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to that at LHC ( 0.007 ^ r ^ 0.07) (see also fig. 12). Because of this, at Tevatron the 
valence qiiark initiated processes dominate while gluon and sea qnark initiated processes 
dominate at the LHC. As the mass of the sneutrino increases, that is x approaches to 
unity, the i^-factor at Tevatron naturally falls off. At LHC, in the higher mass region (~ 1 
TeV), valence quark densities start to dominate and hence it stays almost flat compared 
to Tevatron. We find that the most dominant sub processes are the dd and dg initiated 
processes. The dg sub process which begins at NLO gives a negative contribution both 
at NLO and NNLO. 

We now turn to study the impact of the factorisation scale {/if) and the renormali- 
sation scale (/ir) on the production cross section. The factorisation scale dependence for 
both LHC (left panel) and Tevatron (right panel) are shown in upper panels of fig. 13, 
for rrii) = 300 GeV (LHC), ma = 120 GeV (Tevatron). We have chosen /ir = mc for both 
the LHC and the Tavatron. The factorisation scale is varied between hf — 0.1 mo and 
/iiT = 10 mj>. We find that the factorisation scale dependence decreases in going from LO 
to NLO to NNLO as expected. 

The dependence of the renormahsation scale dependence on the total cross sections 
for the resonant production of sneutrino at the LHC and the Tevatron is shown in the 
lower panels of fig. 13. Note that the LO is already dependent due to the coupling 
^'{(^r)- We have performed this analysis for sneutrino mass m^ = 300 GeV (LHC), 
= 120 GeV (Tevatron). We have set the factorisation scale iif = rno/A (see for 
example, Ref. [53] for more details) and the renormahsation scale is varied in the range 
0.1 < fJ-R/mp < 10. We find significant reduction in the fXR scale dependence when higher 
order QCD corrections are included. It is clear from both the panels of fig. 13 that our 
present NNLO result makes the predictions almost independent of both factorisation and 
renormahsation scales. 

3.2 Charged slepton production 

We now study the numerical impact of our NNLO results on the charged slepton pro- 
duction for both Tevatron and LHC energies. In fig. 14, we have plotted the total cross 
section as function of charged slepton mass. The upper (lower) set of fines corresponds 
to the factorisation scale hf — mi+{0.1mi+). The improvement due to NLO and NNLO 
pieces is evident. As we expect the factors (see figs. 15) are quite similar to the case 
of neutral scalar production. In upper panels of fig. (16), we have shown the production 
cross section as a function of factorisation scale iif for the slepton mass = SOOGeV 
(LHC), ~ 120GeV (Tevatron) and fixed the renormahsation scale at /ir — m^+. In 
the lower panels of fig. (16), we have shown hr variation for the LHC and the Tevatron. 
We have done this for charged slepton mass = 300 GeV (LHC), mf+ ~ 120 GeV 
(Tevatron) and fixed the factorisation scale fip = \ fT^£+ (see Ref. [53]). We again find 
that the scale dependence gets reduced significantly as we include higher order terms in 
the perturbative expansion. Our numerical code can also produce results for l~ produc- 
tion both at the LHC and the Tevatron. At Tevatron, the production rates for is found 
to be same as that of because the contributing fluxes and the partonic coefficient func- 



13 



tions are identical. At the LHC, the is not case because of different fluxes that contribute 
for i+ and t . 

3.3 Soft gluon contributions at order 

The order partonic coefficient function A^^^ resulting from the soft gluons is known 
for the Higgs boson production through bb annihilation (see [51]). The same coefficient 
function can be used here to study the impact of soft gluons for the sneutrino (also for the 
charged slepton) production at «^ level due to identical structure of the interaction terms 
responsible for their production mechanisms. For our numerical study, we have used the 
available three loop (5^ and 73 to evolve and A' respectively. For parton density sets, we 
can only use the available NNLO evolved MSTW 2008 sets. In fig. (17), we have plotted 
the N^LO corrected sneutrino production cross sections against its mass (upper panels) 
and the renormalisation scale (lower panels) on for both LHC as well as Tevatron energies. 
While the order aij soft gluon effects to the production of sneutrinos is indistinguishable 
with respect to exact NNLO contributions, it reduces the uncertainty resulting from the 
choice of renormalisation scale significantly. 



4 Conclusions 

The ongoing program in search of signals of BSM scenarios can be successful only if the 
theory predictions are precise and reliable. In this paper we have attempted to make pre- 
dictions for the production of sneutrinos and sleptons at hadron colliders such as Tevatron 
and the LHC. The resonant production of sneutrinos and sleptons are possible at these 
colliders thanks to i?-parity violating interactions present in the supersymmetric theory 
which is one of the most studied BSM in the literature. Often predictions based on 
subprocess contributions computed at leading order in perturbation theory suffer from 
uncertainties resulting from the arbitrariness in the choice of renormalisation and fac- 
torisation scales. These scales are artifacts of the perturbation theory and hence are 
unphysical. The sensitivity to these scales signals the missing higher order contributions 
that need to be included in order to make the predictions reliable. In other words, con- 
tributions from NLO and NNLO sub processes are expected to reduce these theoretical 
uncertainties. In addition, the potential SM background processes to resonant production 
of sneutrinos and sleptons and their subsequent decays are well under control as they are 
known to NNLO level in QCD. Hence, we have computed all the subprocess contributions 
to the production cross sections for sneutrino and sleptons upto order (i.e., NNLO) 
in perturbativc QCD. We have used dimensional regularisation to regulate UV, soft and 
coUinear divergences. The renormalisation and the factorisation arc done in M S scheme. 
We have used the latest parton density sets provided by MSTW for all our analysis. We 
have demonstrated how the inclusion of NNLO contributions can reduce the scale un- 
certainties over wide range of sneutrino and slepton masses for both Tevatron and LHC 
energies. We also find significant increase in the cross section due to opening up of several 
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new partonic channels beyond the leading order. The increase in the cross section due 
to the inclusion of NNLO contributions compared to that of NLO is found to be 6.6% to 
21% for Tevatron and 3.4% to 4% for LHC and the total NNLO K factor varies from 1.65 
to 1.85 for Tevatron and 1.27 to 1.52 for LHC in the mass range of 100 GeV to 1 TeV. In 
order to estimate the impact of QCD corrections beyond NNLO, we have performed an 
analysis taking into account the dominant soft gluon contributions at N^LO and found 
that they are stable under perturbation. The calculations presented in this paper are 
not particular to supersymmetric theories, but can be applied to any color-singlet scalar 
(pseudoscalar) coupling to a quark anti-quark pair. 
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5 Appendix 



The leading order contribution (fig.(l)) from the subprocess Qi + Qj ^ (f) gives: 



,(0) _ a(o) 



A ^ = 5{1 - x) 



(38) 



To order a^, we need to include ?i + ?j — >■ 5' + 0, + 5' ^ (Ijilj) + one loop 

corrections to qi -\-q, — > 0. The quark anti-quark initiated processes (fig. (2, 3)) give 
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and quark (anti-quark) gluon initiated processes (fig. (4)) give 
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(40) 



In the above results we have presented results to order e that will contribute to NNLO 
restults. 

To order a^, several partonic subprocesses contribute. We present the contributions 
coming from each subprocess below. The quark anti-quark initiated processes with no 
gluon (fig. (5)), one gluon (fig. (6)) and two gluons (fig. (7)) in the final state along with 
(f) constitute a sub class giving order contribution. This class consists of processes 
with two loop corrections to leading order qi + qj — )> (no gluon) along with one loop 
corrections to + g^- — ?■ + (? (one gluon). In addition we have g^ + g^ —>(? + (? + (two 
gluons) in this class. The unrenormalised form factor J^^ upto two loop order (fig. (5)) is 
found to be [51] 
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The contributions from quark antiquark annihilation processes can be split into two parts: 
contributions coming from threshold region called soft plus virtual (S+V) contribution: 
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The functions Snp and Li^ are Nielsen integral and polylogarithm respectively: 
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and for s = 2,3, we have ^2 = tt'^/G,^^ = 1.202056903159 • ••. In addition, quark anti- 
quark pair in the initial and final states along with (f) (A and B in fig. (8)) and quark 
(anti-quark) quark (anti-quark) pair in the final state with also arise at order (C 
and D in fig. (9)). The contributions coming from the later will depend on whether 
we have identical quarks (anti-quarks) or non-identical quarks (anti-quarks). We present 
their contributions below along with various allowed interferences. 

The s-channel processes (A in fig. (8)) with (p emitted from incoming partons give 
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The interference of s-channel processes (A in fig. (8)) with emitted from incoming partons 
and t-channel processes (C and D in fig. (9)) give 
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The s-channel processes (B in fig. (8)) with emitted from outgoing partons give 



(47) 



CpTf 



16 64 16 2. / X / 32 2 

y + yX - 16X - y X (2 + Ll2 (-x) I - yX 



In (x) I - ^ + y a; + 8xH + In^ (x) | ^x^ 



32 



+ ln(l + x)ln(x) i - yx^ 



(48) 



The interference of s-channel processes (B in fig. (8)) with emitted from outgoing partons 
and t-channel processes (C and D in fig. (9)) give 



{2),BC 



'qq 



6 - 4C2 + 32x + 8xC2 - 26x^ + Si,2 {-x) [- 48x^) 

+5i,2 (1 - x) (32x2) + Lis {-x) ( - 8x2) ^ ^-^ ( - 8 + 16x) 

+Li2 (1 - x) ( - 8 + 32x - 24x2) + In (x) { - 8 + 28x + 8x^(2 

+Li2 (-x) (24x2) + Li2 (1 - x) (l6x2) } + ln2 (x) { - 2 + 12x - 12x2} 



21 



+ ln^ (x) {^a;^} + In (1 + x) { - 24x^(2 + Li2 (-x) ( - A8x'^) 



+ In (x) ( - 8 + 16a;) + In^ (x) {20x^) } + In^ (1 + x) In (x) { - 24:X^} 



(49) 



For charged scalar production, only part of t-channel processes (C or D in fig. (9)) interfere 
with s-channcl processes in of (A and B in fig.(8)), that is, either [A^^'^'" , A^^'^'") or 
(Ag)'^"", A^^j''''^) will contribute. 

To order a^, quark (anti-quark) gluon initiated processes also receive contributions 
namely from one loop contributions to qi{q^) + g — t- qi{qj) + 4> and processes involving an 
additional gluon in the final state, qi{^,j) + (?—?■ qi^Mj) + 5' + 0- This contribution is found 
tobeA(J = Ag = A(|g-+A(|g^where 
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-24a;2) + In (1 - x) |l28 - I6C2 - 394a; + 32a;C2 + 272a;2 - 32a;\2 
+Li2 (1 - a;) ( - 4 + 8a; - 104a;^) + In (a;) (56 - 256a; + 216a;^) 
+ In^ (a;) (24 - 48a; + 80x^) | + In^ (1 - a;) | - 66 + 192a; - 138a;^ 



+ In (a;) ( - 42 + 84a; - 132a;2) | + In^ (1 - a;) | y - + 



-X 



+ ln 



rrij 



54 - 8C2 - 172a; + 16a;C2 + 120a;^ - 16a;^C2 



+Li2 (1 - a;) ( - 48a;^) + In (a;) (22 - 128a; + 116a;^) + In^ (a;) (8 - 16a; 
+32a;^) + In (1 - a;) ( - 72 + 200a; - 140a;^) + In (1 - a;) In (a;) ( - 40 

+80a; - 128a;^) + In^ (1 - a;) (36 - 72a; + 72x'^) | + In^ | - 15 

+36a; - 24a;^ + In (a;) ( - 6 + 12a; - 24a;^) + In (1 - x) 



X (12 - 24a; + 243;^) | 



(51) 



We now present the contributions resulting from pure i-channel processes where the final 
state quarks (anti-quarks) are non-identical (C and D in fig. (9)). The individual i-channel 
results are found to be identical and are given by 



A {2),CC ^ a{2),DD ^ A (2),CC ^ a{2),DD ^ A (2).CC ^ a (2),DD 
qq qq qq qq qq qq 



C^Tf 



305 208 32 



470 



- ^C2 - 8C2 - + 8XC2 + + —X\2 



703 2 32 

x^ H 

27 3 



9 27a; 3a; 

+5i,2 (1 - x) (16 + 16a;) + Li3 (1 - a;) ( - 32 - 32a;) + Lis (1 - x) 

176 



/ 64 32 .\ , , 62 

x( - 4+ — + 28a; + ya;M + In (a;) y - I6C2 - — 



-X — 16x^2 



+ + Lia (1 - a;) ( - 8 - 8a;) I + In^ (a;) | - 1 - 5a; - 



24 



, 3 / ^ I 10 10 1 , ,^ , [ 152 32 200 176 ^ 

+ In-^ {x){ \ X ) + hi{l- x) { \ \ X x^ 

^ '\ 2> 31 ^ '\ 3 9a; 3 9 



+Li2 (1 - x) (32 + 32a;) + In (a;) ( - 8 + 32a; + 32a;^ 



32 


32 2 






3a; 


8a; — ^a; 


16 


100 88 2 




X X 


to ^ 


3 9 



+ ln(x) (l6+16x)|+ln(^)|-| 
+Li2 (1 - a;) (l6 + 16a;) + In (a;) ( - 4 + 16a; + 16a; 

(oo oo \ 

8 + — - 8a; - —x"^ \ 

+ In (1 - a;) In (a;) (l6 + 16x) | + In^ j |2 



+ In (a;) (4 + 4a;) 



(52) 



If the final state quarks (anti-quarks) are identical, we find in addition to above men- 
tioned t-channcl processes, we have to include similar t processes with final state quarks 
(anti-quarks) interchanged (fig.(lO)). Hence we have contributions from the additional t 
processes (fig. (10)) and their interferences with the already existing t-channel processes 
(fig.(9)). The former gives results identical to A^^)'*^'-^, A^^^'^^ and A^^)'*-^^. The interfer- 
ence contributions are found to be 



A (2),CE ^ A (2),DF ^ A (2),CiJ ^ A 
qq qq qq qq 



1 — a; 
1 + x 



-34- 



l + x 



Cs + 4(3 + 4(2 + 34a; - 4a;(3 



- 32 + + 32a;^ + Lis f-^) { 32 

1 + a; ) \ l + xj \ 

^ - 32a;j + S^,^ {-x) (^40 - ^ - 40a;j + 5i,2 (1 - x) 
X I - 32 + + 32a;) + Li^ {-x) (12 - -r^ 12a; 



l + x 



l + x 



+Li3 (1 - a;) 1^32 - - 32a;j + Li2 {-x) (s) + Li2 (1 - a;) 



25 



X 



{24 
- 18 + Y^C2 - I6C2 + 14a; + 16a;C2 



64 48 
+Li2 i-x) (-44 + + Ux) + Li2 (1 - x) (-24 + + 24x) 



+ In^ (x) { - 4 - lOx} + W (a;) I y - 



10 16 10 1 , 

o"-^ > + In (1 + 



31 + a; 3 



1 16 / 32 ^ 

X <^ - C2 + 2OC2 - 2OXC2 + Li2 i-x) 40 - 40x 

\ 1+X \ 1+X J 



+ In {x) (8) + In^ (x) f - 38 + + 38x j I + In^ (1 + x) In {x 



( 16 1 (32 

X <^ 20 20a; ^ + In (1 - x) <^ 32 (2 + I6C2 - 32a; 

1 + a; 1 + a; 



-16a;C2 + Li2 (-x) i 32 - 32a; j + In (a;) (16 + 16a;) 



1 + a; J 



+ In^ (a;) ( - 8 + + 8a; ) + In (1 + x) In (a;) (32 - - 32a; 



1 + X 

16 



1 + a; 



+ In ( -fj |16 ^ + 8C2 - 16a; - 8a;C2 + Li2 {-x) I 16 



32 



1 + a; 



- 16a; j + In (a;) (8 + 8a;) + In^ (a;) ^ - 4 + 



1 + a; 



+ 4a; 



+ In (1 + a;) In (a;) fl6- - 16a; 



A (2),Dfi ^ A (2),CF ^ A (2),DB 

QQ m m 



23 + 36a; - 13a;2 + 81^2 (1 - a;) (4 - 8a; - 8a;^) 

+Li3 (1 - a;) (4 - 8a; + 8a;^) + Li2 (1 - a;) (4 + 12a; - 12a;^) 

+ ln (a;) {2 - 4a; + Li2 (1 - a;) ( - 8a;^)} + In^ (a;) (2 + 6a; - 6a;^) 



26 



(55) 



Finally wc compute gluon gluon initiated processes that contribute to order a^. Writing 

^(2) ^ ^{2),Ca + /^{ihCj. 
99 99 99 ' 

N"^ \ 1 , 248 
3+8(3-^ 



A(2),Ca ^ 

99 A/"2 - 1 



+^1,2 (-a;) (8 + 16a; + 16x^) + ^1,2 (1 - x) ( - 4 + 8a; - 8a;2) 
+Li3 (-x) (12 + 24a; + 24a;^) + Lia (-a;) ^ - 8a; + ya;^ j 

+ In (x) 1^ - ya; - 58a;^ + Lis (-x) ( - 12 - 24a; - 24a;^) | 

+ In^ (x) |2x + yx^l + In (1 + x) {4(2 + 8^(2 + 8x^(2 

+Li2 (-x) (S + 16x + 16x^) + In (x) ^ - 8x + yx^ j + In^ (x) 

+ In^ (1 + x) In (x) |4 + 8x + 8x^| 



X - 6 - 12x - 12x^ 



(56) 



and 



A(2),Cj. 
99 



Si,2 i-x) ( - 8 - 16x - 16x^) + Si,2 (1 - x) ( - 12 - 72x - 56x^) 
+Li3 (-x) ( - 12 - 24x + 8x2) ^ j^-^ _ j^^g ^ ^ g4^2^ 
+Li2 (-x) (8x) + Li2 (1 - x) ( - 4 + 16x + 56x2) _^ ^^-j 
X I - 15 + 8C2 - 48x + 32xC2 + 121x2 + 40x^^2 + Li2 (-x) (l2 

+24x + 8x2) ^ 2>i2 (1 - x) ( - 4 - 16x - 16x2) | + ln2 (x) { - 4 

{2 8 16 1 f 
- - - -X - yx2 I + In (1 + X) <^ - 4(2 

-8XC2 - 8x2^2 + Li2 (-x) ( - 8 - 16x - 16x2) + In (x) (8x) 

+ ln2 (a;) (6 + 12x + 12x2) | + ln2 (1 + x) In (x) { - 4 - 8x - 8x2} 
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+ In (1 - x) {46 + lOAx - 150x^ + Li2 (1 - x) ( - 16 - Ux - 64x^) 

+ In {x) (20 + 64x - 16x2) + In^ {x) (a + 16x + IQx"^) } 

+ In^ (1 - x) { - 16 - 32x + 48x2 + In (x) ( - 8 - 32x - 32x2) } 

+ In j |23 + 52x - TSx^ + Lia (1 - x) ( - 8 - 32x - 32x2) 
+ In (x) (10 + 32x - 8x2) ^ ^^2 ^^-^ (2 + 8x + 8x2) + in (1 _ 
X ( - 16 - 32x + 48x2) + In (1 - x) In (x) ( - 8 - 32x - 32x2) | 

+ ln2 j I _ 4 _ 8a; + 12x2 + In (x) ( - 2 - 8x - 8x2) | 

-20 - 8C3 + I6C2 - 98x - I6XC3 + 36xC2 + 118x2 + 8x2(3 - 48x2(2 (57) 
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Figure 1: Subprocess + qj — > (f). 
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Figure 3: Subprocess qi + qj ^ (j) + g. 
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Figure 4: Subprocess qi + g ^ (j) + qj . 
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Figure 5: Subprocess qi + qj — )■ 0. 
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Figure 6: Subprocess + — )■ + g. 
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Figure 7: Subprocess qi + qj ^ (j) + g + g. 
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Figure 9: t channel gluon exchange graphs contributing to the subprocesses qi + qj 
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Figure 10: t channel gluon exchange graphs contributing to the subprocesses q{q) + q{q) 
(p + q{q) + q{q) with identical quarks (anti-quarks). 
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Figure 11: Total cross-section for the v 'production. The upper (lower) set of lines corre- 
spond to the factorisation scale fip = m^(0.1mj>). 
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Figure 12: K-factor for the v production. 
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Figure 13: fip and /i^ variations for the v 'production. 
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Figure 14: Total cross-section for the production. The upper (lower) set of lines 
correspond to the factorisation scale fip = mi+{0.1mj^) . 
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Figure 15: K-j actor j or the 'production. 
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Figure 16: fip and variations for the £^ production. 
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Figure 17: Total cross section and fin variations for the v 'production. 
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